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Abstract: Although little morphological differentiation is detected between intercontinental disjunct populations of Phryma, 
this monotypic genus shows distinct molecular divergence corresponding to its distribution in eastern Asia and eastern North 
America . This study further employs amplified fragment length polymorphism (AFLP) analyses using two selective primer 
pairs to quantify genetic diversity and evaluate the morphological stasis . Most of the molecular variance is accounted for by 
variance among populations between regions . Cluster and PCA analyses revealed that Phryma constitutes two major groups in line 
with their geographic distribution, with one genetically distinct group from the eastem Asia and the other with accessions from 
eastern North America . The results robustly supported the distinct genetic divergence and morphological stasis in Phryma . 
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Phryma L . ( Phrymaceae) is well - known for its opmentally) and the classical intercontinental disjunct 
unusual morphology with a derived pseudomonomerous distribution between eastern Asia and eastern North 
gynoecium (bicarpellate with one carpel reduced devel- America (Li, 1952; Wen, 1999) . It is also one of the 
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very few monospecific genera with such an interconti- 
( Hara, 1966; 
Thieret, 1972) or subspecies (Kitamura and Miurata, 
1957; Li, 2000) have been recognized within the sin- 


nental disjunction. Two varieties 


gle intercontinental disjunct species Phryma lep- 
tostachya L ., corresponding to its biogeographic distri- 
bution in eastern Asia and eastern North America . A 
high level of morphological similarity has been reported 
for the two disjunct subspecies ( Hara, 1969; Ramana 
et al., 1983), but companied with distinct molecular 
divergence (Nie et al., 2006) . 

Hara (1966) pointed out that the intercontinental 
disjunct populations of Phryma were identical in most 
cytology, 
habitats . For example, populations from both conti- 


morphological characters, and ecological 
nents bear flowers that are erect in bud but later 
spreading or becoming deflexed; they have the same 
chromosome number of 2n = 28; and they share similar 
habitats of deciduous or mixed forests . Plants from the 
two different regions only differ slightly in leaf size, 
shape of upper lip of the corolla, and length of the up- 
per spinulose calyx-lobes (Hara, 1962, 1966, 1969; 
Li, 2000; Nie et al., 2006) . 

Although a high level of morphological similarity 
has been found in Phryma, recent molecular studies 
revealed substantial molecular divergence between the 
two intercontinental varieties ( Beardsley and Olmstead , 
2002; Lee et al., 1996; Nie et al., 2006; Xiang et 
al., 2000) . A low level of isozyme identity (0.291) 
from allozyme loci and high ITS divergence (4.46% ) 
were detected between the disjunct pair in Phryma 
(Lee et al., 1996) . 

The eastern Asian-eastem North American dis- 
junction is a well-known and classical biogeographic 


pattern in theNorthern Hemisphere that has received 


Table 1 
Taxon Abbreviation 
EN1 
Phryma leptostachya L . var. leptostachya EN2 
(eastern North America) EN3 
EN4 
Phryma leptostachya var. asiatica Hara EAI 
(eastern Asia) EA2 
EA5 


considerable attention in the last 15 years (Donoghue et 
al., 2001; Meng et al., 2003; Milne, 2006; Milne and 
Abbott, 2002; Ran et al., 2006; Wen, 1998, 1999, 
2001; Wen et al., 1996; Wen et al., 1998; Wen and 
Stuessy, 1993; Wen and Zimmer, 1996; Xiang ef al., 
1996, 1998, 2000; Zhou et al., 2006) . Most studies 
on taxa with disjunct distribution in these two regions 
(especially for the disjunct species) have focused on 
phylogenetic relationship and biogeographic history at 
the interspecific level (Nie et al., 2006), yet few 
have analyzed patterns of genetic variation among the 
disjunct taxa at the population level . The objectives of 
this preliminary study are to explore genetic divergence 
between intercontinental disjunct populations of Phry- 
ma and examine the phylogeographic structure and bio- 
geographic implications of Phryma in the Northern 


Hemisphere . 


1 Materials and Methods 

Seven populations were analyzed with three from eastern 
Asia and four from eastern North America (Table 1) . Genomic 
DNA was extracted from 15 mg silica dried leaf material using the 
modified CTAB method (Doyle and Doyle, 1987) . 

Procedures for AFLPs were as described by Vos et al. with 
some modifications (Vos et al., 1995) . 400 - 500 ng DNA of 
each sample was digested and ligated with 5 U of EcoRI and 1 U 
of Msel in a reaction mixture that contained 1X T4 restriction and 
ligase buffer with ATP, 0.5 Umol L`’ EcoRI adapter, 5 pmol 
L`! Msa adapter, 50 mmol L` ! NaCl, 5 ng BSA, 1 U T4-DNA 
ligase, and water to a final volume of 10UI1 . After digestion and 
ligation, products were diluted 1: 20 with distilled water . In the 
first round of PCR (pre-selective amplification), each reaction 
contained 2.0 ul Promega 10X reaction buffer, 2.0U1 dNIPs (2 
mM), 10.0 Ul water, 2.0 ul Promega Taq Polymerase (5 U 
UL’), 0.51 Msel+C primer (50 nguL'), 0.51 EcoRI 
+ A primer (50 nguL`'), and 2.0u1 of the diluted digestion 


Voucher information of Phryma populations sampled in this study 


Voucher Locality 
Wen 7140 (US) USA: Illinois 
Wen 7161 (US) USA: Illinois 
Wen 7188 (US) USA: Alabama 
Wen 7292 (US) USA: Wisconsin 


Nie 102 (KUN) China: Yunnan 
Yue 122 (KUN) 


Wen 5757 (US) 


China: Sichuan 


China: Yunnan 
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ligation product . The PCR conditions were: 94°C for 20 sec, 
56°C for 30 sec, and 72°C for 2 min, for 19 cycles . The PCR 
products were diluted 1: 20 with water . Samples used in the 
bulked DNA study were run individually through the preselective 
amplification step . Subsequently, 10 ul of the diluted preselec- 
tive amplification from each of three individuals was combined . 
In the second round of PCR (selective amplification) , each reac- 
tion contained 1.51 Promega 10X reaction buffer, 1.0u1 dNTPs 
(2mM), 11.2751 water, 1.51 Promega Taq Polymerase (5 U 
UL’'), 0.5 Ul Msel + XXX primer (50 nguL’'), 0.5 ul 
EcoRI + XXX (+ dye) primer (50 nguL™'), and 1.51 of the 
preselective product . The PCR conditions were: 94°C for 20 
sec, 66°C for 30 sec, and 72°C for 2 min, 9 cycles; 94°C for 
20 sec, 56°C for 30 sec, and 72°C for 2 min, 30 cycles . All 
samples were processed in random order and subsamples were re- 
run for internal control . Selective amplification products were 
separated using Long Ranger Singel packs ( BioWhittaker Molecu- 
lar Applications, cat . # 50691) running for about 2 hours on 
the ABI Prism 377 genetic analyzer by combining 0.4 Ul of selec- 
tive amplification product with 1.2u1 of loading buffer containing 
0.751 formamide, 0.30 ul GeneScan ROX-500 internal size 
standard, and 0.15Ul loading dye . 

We identified and sized the peaks between 50 and 500 bp in 
the ABI gel images using GeneScan 3.1 (Applied Biosystems ) to 
create sample files . Sample trace files were imported into Genog- 
rapher (ver 1.6, Montana State University, 2001 ; http -@horde- 
um.msu .montana. edæenographe® for visualization and scor- 
ing . AFLP loci were analyzed using the“ thumbnail” option of 
Genographer, which allows for comparison of signal strength at 
each locus for all samples . Only polymorphic AFLP markers were 
scored for each sample and recorded as a binary character with 


LL 


presence as“ 1” and absence as“ 0” . Monomorphic markers 
were not scored . 

The genetic differentiation between populations of Phryma 
based on all AFLP loci was calculated using Nei s genetic dis- 
tance (Nei, 1978) . A cluster analysis using the unweighted pair- 
group method with arithmetic averaging ( UPGMA ) (Sneath and 
Sokal, 1973) was performed using the software POPGENE 1.32 
(Yeh etal., 1997) .A genetic distances matrix was also used to 
perform a hierarchical analysis of molecular variance (AMOVA) 
(Excoffier et al., 1992), using ARLEQUIN version 3.11 ( Ex- 
coffier et al., 2005), to partition the total genetic variation 
among individuals within populations, between populations within 
a region, and between regions of eastern Asia and eastern North 
America . These are not directly comparable to the traditional F- 
statistics (Wright, 1951) and are therefore referred to as O-sta- 
tistics (Excoffier ea al., 1992) . The significance of Ọ values 
was tested by 1000 permutations . 


The AFLP data were also subjected to a principal compo- 
nents analysis (PCA) , which may help reveal unexpected rela- 
tionships among a large number of variables into two or three new 
uncorrelated variables so that they retain most of the original in- 
formation . Similarity matrices using the Jaccard s similarity coef- 
ficient were generated . Eigenvalue and eigenvector matrices were 
calculated from the similarity matrix . The standardized data were 
projected onto the eigenvectors of the correlation matrix and repr- 
esented in a two-dimensional scatter plot . Plots of samples in re- 
lation to the first three principal components were constructed 
with populations designated as either eastern Asian or eastern 
North American distribution . PCAs were performed using the 
computer program NTSYSpc, version 2.11h (Rohlf, 2000) . 


2 Results 

The two AFLP primer pairs M-CT@@E-ACT and 
M-CT@E-AGG generated a total of 426 polymorphic 
markers from 80 individuals for seven populations . Ta- 
ble 2 shows the genetic differentiation [Nei genetic 
distance; (Nei, 1978)] between Phryma populations . 
The genetic distance matrix was used to establish the 
level of genetic divergence between the populations 
(Table 2) . Estimates of genetic distance using AFLP 
data ranged from 0.027 for the most closely related 
populations within each region ( EN1 and EN2), to 
0.811 in the most divergent populations between re- 
gions (EAI and EN3) . These results indicate that ge- 
netic distance is solely dependent on geographical dis- 
tance, and clustering based on genetic distances re- 
flects geographical relationships . 

Table 2 Nei s genetic distance of Phryma populations estimated 
by AFLP analysis (Population codes are as in Table 1) 


Population EAI EA2 EAS ENI EN2 EN3 EN 


EA2 0.1367 - 

EA5 0.1510 0.1248 - 

EN1 0.6233 0.5133 0.4235 

EN2 0.6214 0.5438 0.4379 0.0274 

EN3 0.8114 0.6483 0.5826 0.0496 0.1032 

EN4 0.6526 0.5948 0.5361 0.0689 0.0535 0.1179 


The results of the AMOVA based on AFLP varia- 
tion are shown in Table 3 . Most of the total variance 
was found between regions (60.57%, P<0Q.0001) . 
The corresponding ®cr value was 0.61 ( P<0.0001) . 
The among-population@wvithin-regions percentage was 
relatively low (18.2% , P<0.0001) when compared 
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Table 3 Results of hierarchical analysis of molecular variance (AMOVA) based on AFLP markers between eastern North American and eastern 


Asian populations of Phryma . P-values, calculated from a random permutation test (1000 replicates), and Ọ -statistics represent the 


probability of obtaining by chance alone a more extreme variance than the observed values ( Excoffier et al., 1992) 


Source of variation d.f. SS Variance component % of variance Ọ-statistics P 
Between regions 1 396.084 9.38894 60.57 Dcr = 0.60574 < 0.0001 
Among populations within regions 7 187.236 2.82036 18.20 Doc = 0.46152 < 0.0001 
Within populations 73 240.216 3.29064 21.23 ® sr = 0.78770 < 0.0001 
Total 81 823.537 15.49993 


to the within-population variance component, close to 
that of the among-individuals 
(21.23%, P <0.001) . All three Ọ-values were sig- 
nificant based on 1000 permutation tests (Table 3) . 
The UPGMA dendrogram based on the genetic 


distances between populations showed two main clusters 


within populations 


of populations that correlated to their continental distri- 
bution (Fig .1) . The PCA of the AFLP-based distance 
data was performed to examine relationships among 
Phryma populations . The first and second principal co- 
ordinates described approximately 44% and 11% of 
the total variation, respectively (Fig .2) . Similar to 
the UPGMA dendrogram, the PCA showed that eastern 
Asian populations clustered apart from those of eastern 
North America . 


3 Discussion 

The AMOVA analyses showed that the largest por- 
tion (61% ) of genetic variance is contributed by ge- 
netic variation between intercontinental regions of east- 
em Asia and eastern North America and only a small 
part is due to divergence among populations within re- 
gions (18% ) and within populations (21%) . The re- 
sults are similar to the pairwise sequence divergence of 
ITS, rpsl6, and trnL-F markers (Nie et al., 2006) . 
For example, the ITS divergence of the two varieties 
between eastern Asia and eastern North America 
(3.11%-4.41%) was higher than that among popula- 
tions within each continent (1.63% - 0.65% ) . The 
high level of genetic divergence between regions sug- 
gests that variations between the two continents display 
a large proportion of genetic loci of high allelic varia- 
tions . 

The long time of geographic isolation (Lee et al., 
1996; Nie et al., 2006; Xiang et al., 2000) may 
contribute to the high molecular differentiation of Phry- 


ma between regions observed . Several previous studies 
estimated the divergence times for the disjunct Phryma 
varieties using various dating approaches . Lee ef al. 
(1996) reported the divergence time to be over 20 mil- 
lion years ago (mya) using allozyme data and 12.35 mya 
using ITS sequences . Xiang et al. (2000) estimated the 
divergence times for the genus with rbcL sequence data 
and a molecular clock calibrated with Cornus fossils . 
They estimated that the two varieties of Phryma di- 
verged about 5.85+ 2.66 mya . Recently we estimated 
their divergence time as 3.68+ 2.25 to 5.234 1.37 mya 
based on combined chloroplast data using both the 
Bayesian dating and the penalized likelihood methods 
with relaxed molecular clocks (Nie et al., 2006) . 
Genetic divergence within the eastern North 
American populations is less than that within the east- 
em Asian populations (Table 2), although only a few 
populations sampled in this study . A similar situation 
is found with the DNA sequences with more extensive 
sampling from our previous study (Nie et al., 2006) . 
The eastern Asian populations are reported to be more 
heterogeneous in ITS sequence variation than those in 
the eastern North America (maximum divergence of 
1.63% vs . 0.65%) . As suggested by Nie et al. 
(2006), the more genetic variation in the eastern 
Asian Phryma may be due to higher level of geographic 
isolation with the more heterogeneous and pronounced 
topographies in easter Asia ( Qian and Ricklefs, 
2000; Wen, 1999; Xiang et al., 2004) . 
Morphological stasis in Phryma is also supported 
by genetic analysis at the population level . Little mor- 
phological variation has been observed from the disjunct 
pairs of Phryma (Holm, 1913; Li, 2000) . Multivari- 
ate analysis based on 23 quantitative morphological 
characters also revealed that no significant variations 
were found between them (Nie ef al., 2006) . However, 
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Fig .1 UPGMA dendrogram based on pairwise genetic distance showing relationships between Phryma 
populations from AFLP data . (@ = EN samples; © = EA samples) 
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Fig .2 Principal component analysis of AFLP data . Accessions are 
plotted according to the values of the first (x-axis) and the second 
(y-axis) components and with different symbols according to 


geographical origin (@ =EN samples; © =EA samples) 


AFLP analysis revealed a clear geographical pattern of 
genetic variation, consistent with our findings from the 
sequence data (Nie et al., 2006) . The results showed 
that Phryma is a genetically diverse entity, and AFLP 
markers can be effectively employed to assess genetic 
diversity and to measure genetic relationship among 
populations . 

Morphological stasis was suggested to explain the 
discordance of the genetic and morphological rates of ev- 
olution (Nie et al., 2006) . Among various possible ex- 
planations for stasis in morphology, a relatively constant 
environment with the concomitant action of stabilizing 
selection, might be the most plausible (Nie ef al., 
2006; Wen, 1999) . The intercontinental populations of 
Phryma occupy similar habitats in rich mesic to moist, 
deciduous or mixed deciduous and evergreen forests in 
both eastern Asia and eastern North America . These 
types of forests covered most of the temperate regions of 
the Northern Hemisphere in the Tertiary, but can be 
found today only in southeast North America and in 
eastern to central China and central to southern Japan 
(Milne, 2006; Milne and Abbott, 2002; Nie et al., 
2008; Wen, 1999) . 

In present study, our sampling is limited with 
only three eastern Asian populations from SW China 
and four populations from eastern North America . Nev- 
ertheless, our AFLP results at the population level cor- 


roborate our findings in the sequence analyses (Nie et 
al., 2006) . Both the AFLP (Table 2) and DNA se- 


quence data (Nie et al., 2006) suggested that the 
eastern Asian populations have much higher genetic and 
molecular divergence than that of eastern North Ameri- 
( Avise, 


2000) of the disjunct taxa is very important to our un- 


ca. Comparative phylogeographic analysis 


derstanding of mechanisms responsible for their phylo- 
genetic relationships and disjunct distribution . Future 
studies are planned to extensively sample the popula- 
tions on both continents to better understand the rela- 
tionships among genetic divergence, morphological dif- 
ferentiation, and the development of the intercontinen- 
tal disjunct pattern . We intend to employ phylogeo- 
graphic analysis, as implemented in Gao et al. (2007) 
and Meng et al. (2007), but analyze the Phryma 


populations intercontinentally and intracontinentally . 
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